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Abstract 

We derive here the Friedland-Tverberg inequality for positive hyperbolic polyno- 
mials. This inequality is applied to give lower bounds for the number of matchings in 
r-regular bipartite graphs. It is shown that some of these bounds are asymptotically 
sharp. We improve the known lower bound for the three dimensional monomer-dimer 
entropy. We present Ryser-like formulas for computations of matchings in bipartite 
and general graphs. Additional algorithmic applications are given. 
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1 Introduction 

The aim of this paper is to explore the connections between the problem of counting the 
number of partial matchings in graphs and positive hyperbolic polynomials, and Ryser-like 
formulas for partial matchings. Given a graph G = (V, E) on N vertices, i.e. = N, 
we want to compute the number of m-matching, i.e. the number of subsets M of edges E, 
where #M = to, and no two edges in M have a common vertex. 

Our main results are for bipartite graphs G := (V\ U \<2,E), where E C V\ x V2 and 
n = #Vi = #V2. Let A(G) G {0, l} nx ™ be the incidence matrix of the bipartite graph 
G. Then the number of 771-matchings in G is equal to perm m A(G), where perm m A is the 
sum of to x m minors of A G K nx ™. For m = n, perm A(G), the permanent of A{G) is the 
number of perfect matchings in G. 

It is well known that the computation of the number of perfect matching in a general 
bipartite graph is j^P — complete. See [32] for the first proof and [3] for a simplified proof. 

Ryser's algorithm to compute the permanent of any A G R raxn [28] remains the most 
efficient exact algorithm, even though it uses around n2 ra operations. One can speed up 
significantly the approximate computation of perm A. One knows to compute the permanent 
of a nonnegative matrix A = [a^] G M™ x " within a simply exponential factor [2]. In the 
case that all the entries of the matrix are uniformly bounded below and above by positive 
constants, one can improve the estimates of the exponential errors [26] and [15]. A fully 
randomized polynomial approximation scheme frpas for the number of perfect matchings 
in a bipartite G, and more generally for the permanent of a nonnegative matrix is given in 
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[25]. This result is generalized in [12] to the number of m-matchings in bipartite G, and 
more generally to perm m A, for any A £ n . 

We now describe our main results for pcrm m A, where A is doubly stochastic, and their 
applications to lower bounds on partial matchings in bipartite graphs. Recall that the 
minimum of the permanent of n x n doubly stochastic matrices, denoted by f2„, is achieved 
only for the flat matrix J n , whose all entries equal to — . Thus perm B > perm J n — for 
any B G Sl n and this inequality was conjectured by van der Waerden [33]. This conjecture 
was independently proved by Egorichev [5] and Falikman [7]. We call the above inequality 
Egorichcv-Falikman-van der Waerden (EFW) inequality. The asymptotic behavior of EFW 
incquaity is captured by the inequality perm B > e~" for any B £ Q n . This inequality 
was shown by the first name author [8] three years before [5, 7]. Let F(n, r) be the set of 
all r-rcgular bipartite graphs G on 2n vertices. For G £ r(?i, r) the matrix B := ^A(G) is 
doubly stochastic. Hence the number of perfect matchings in G is at least (-)". Thus for 
r > 3, the number of perfect matchings in r-regular bipartite graphs grows exponentially, 
which proves a conjecture by Erdos-Renyi [6]. Schrijvcr [29] improved EFW inequality to 
r-regular bipartite graphs, whose asymptotic growth is best possible. Recently, the second 
name author [19] improved Schrijvcr's inequality. Moreover, the proof in [19] is significantly 
simpler and transparent. One of the main tools in the proof in [19] is the use of the classical 
theory of hyperbolic polynomials. 

It was shown by the first name author that perm m A > perm m J n for any A £ 
and for m £ [2,n] equality holds only if and only if A = J n [9]. (perm X A = n for each 
A £ Q n .) This was fact was conjectured by Tverberg [31], and is called in this paper the 
Friedland-Tverberg (FT) inequality. FT inequality gives a lower bound on the number of 
partial matchings in any G £ T(n,r). 

We derive here the Schrijver type inequalities for m matchings in r-regular bipartite 
graphs on 2n vertices. This is done using the results and techniques of [19]. In particular we 
give a generalized versions of FT inequality to positive homogeneous hyperbolic polynomials, 
which are of independent interest. 

The notion of partial matching in r(n,r) can be extended to asymptotic matchings as 
n — > oo as follows. Given a sequence of Gi £ r(nj, r) we can consider the quantities 

tt m\\ r . logPerm m ,4(G|) log perm A(Gi) 
f(p,{Gi}) :=hmmf — '■ , F(p, {Gi}) := hmsup — '■ , (1.1) 

m->oo Zni n,->-oo 2ni 

TThl 

where ni — > oo and lim — = p £ [0, 1]. 

f(j>, {Gi}) and F(p, {G{\) can be viewed as the minimal and the maximal exponential growth 
of matchings of density p of the sequence Gi,l £ N. 

Consider the following special case of the above example. Let C m be a cycle on m 
vertices. Note that C m is bipartite if and only if m s even. Fix a positive integer d 
and let Tii t d '■= C21 x ... x C21, be bipartite toroidal grid on (2l) d vertices. Note that 

d 

T 2 i,d e T{2 d - 1 l d 1 2d). It is shown in [21] that f{p,{T 2l>d }) = F(p,{T 2lA }) and this quantity 
is the exponential growth rate of the number of monomer-dimcr tilings of the d-dimcnsional 
cubic grid having sides of length as n tends to infinity and the dimer density, (fraction of 
the maximum possible number ^n d of dimers), in these tilings converges to a fixed number 
p £ [0, 1]. (See also [14].) We denote by this exponential growth by hd{p), and call it the 
c?-dimcnsional monomer-dimer entropy of dimer density p £ [0, 1] in the lattice Z d , see [21] 
and [13]. For d = 1 the rate is known explicitly as a function p, whereas for d > 1 the exact 
rate is unknown. 

hd(p) can be estimated if one can estimate the quantities /(p, {Gi}) < F(p, {Gi}) from 
below and above for any sequence Gi £ Y(ni,r),l £ N. Lower and upper estimates of 
f(p,{Gi}) < F(p,{Gi}) are conjectured in [10] and called the asymptotic lower match- 
ing conjecture and asymptotic upper matching conjecture, abbreviated here by ALMC and 
AUMC respectively. For r = 2 ALMC and AUMC are proved in [11]. 
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In this paper we apply our lower bounds on the sum of subpermanents of doubly stochas- 
tic matrices with r nonzero elements in each row to obtain lower bound on f(p, {G;}). For 
a fixed integer r > 3, we show the validity of ALMC for the densities p s = where 
s = 0, 1, . . .. These inequalities yield new lower bounds for the d-dimensional monomer- 
dimer entropy of dimer density h<i(p),p £ [0, 1] in the lattice Z d . In particular we obtain 
the best known lower bound for the three dimensional monomer dimer entropy which 
combined with the known upper bound in [13] gives the tight result h 3 £ [.7845, .7863]. 

Next we discuss briefly the sum of m x m subhafnians of 2n x In symmetric B with 
nonncgative entries, denoted by haf m B. For 0—1 matrix B this is equivalent to the 
number of m-matchings in a general graph on 2n vertices. Wc give Ryscr-type algorithm 
for perm m A and haf m B. Unfortunately the generating function x T Bx, (the quadratic 
for associated with B), is positive hyperbolic if and only if the second eigenvalue of B is 
nonpositive. We show that for any graph G, x T £?(G)x is positive hyperbolic if and only 
if G is a complete fc-partite graph. The last section is devoted to algorithmic applications 
related to perm m A and haf m B. 

We now list briefly the contents of this paper. In §2 we discuss briefly the notion of 
positive hyperbolic polynomials and their properties that needed here. In §3 we bring the 
generalized version of FT inequality for positive hyperbolic polynomials. In §4 we state and 
discuss the ALMC and AUMC. In §5 we give lower bounds on /(p, {G;}). We apply these 
bounds to verify the ALMC for a countable values of densities for each r > 3 as explained 
above. In §6 wc discuss the notion of haf m B and its connection to the quadratic form 
x T Bx. In §7 we discuss the Ryser-type formulas for perm m A and haf m B. 

We thank Uri Peled for supplying us with the Figures 1 and 2. 

2 Positive hyperbolic polynomials 

Definitions and Notations 

1. A vector x := (xi, . . . ,x n ) T £ R™ is called positive or nonnegative, and denoted by 
x > := (0, . . . , 0) T or x > if Xi > or X{ > for i = 1, . . . , n respectively. A 
nonnegative vector x ^ is denoted by x ^ 0. y > x y — x > 0. The cone of 
all nonnegative vectors in R™ is denoted by R™ . 

2. A polynomial p = p(x) = p(x\, . . . , x n ) : R™ — > R is called positive hyperbolic if the 
following conditions hold: 

• p is a homogeneous polynomial of degree m > 0. 

• p(x) > for all x > 0. 

• cj)(t) := p(x + tu), for iGR, has m-real i-roots for each u > and each x. 

3. For any polynomial p : I" — ► R and any 0/u = (u\, . . . , u n ) T £ R™ let p u = p u (x) : 

4. Let e; := (8n, . . . , Sin) 1 £ R™, i = 1, . . . , n be the standard basis in R". 

5. Let 1 :— (1, . . . , 1) T £ R" and denote by J n £ R nx " the n x n matrix whose all entries 
are equal to i. 

The following lemma summarizes the basic properties of positive hyperbolic polynomials 
that needed here. 

Lemma 2.1 Let p : R™ — > R be a positive hyperbolic polynomial of degree m > 1. Then 
the following properties hold: 
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1. Let u ^ 0,x be fixed and denote <p(t) = p(x + tu). Assume that p(u) > 0. Then <fi(t) 
has m real t roots. Furthermore p u (x) is a positive hyperbolic polynomial of degree 
m—1. y > x > ^ p(y) > p(x) > 0. 

2. Let u ^ 0,x ^ and assume that p(u) = 0. Then either (f>(t) > for all t > or 
p(x) = and 4>{t) = 0. Assume that p(x) > and <p{t) is not a constant polynomial. 
Then all its roots are real and negative. If p u is not a zero polynomial then p u is a 
positive hyperbolic of degree m—1. 

3. If q{{x\, . . . , x„_i)) := p((x±, Xi, ■ ■ ■ , ^n-lj 0)) is not identically zero then q is a posi- 
tive hyperbolic of degree m in R™~ . In particular, r((x±, . . . ,x„_i)) := 

-g2-{(xi, . . . , X n -\, 0)) is either a zero polynomial or a positive hyperbolic polynomial 
in n — 1 variables of degree m — 1. 

Proof. 

We assume that n > 1 otherwise all the results are trivial 

1. The facts 4>(t) has m negative roots and p u is positive hyperbolic is in [17] or [24]. 
Hence p u > on R™. Let v > and assume that u; > 0, 1 = 1, . . . be a sequence 
of vectors converging to v. Then p Ul — > p v . Therefore p v > on R™. In particular 
p ei (x) = J!f-(x) > n R™ for i = 1, . . . , n. Thus p(x) is a nondecreasing function in 
each variable Xi on i?" . Hence y > x > p(y) > p(x) > 0. 

2. Recall Brunn-Minkowski theorem that / := p™ is convex on all positive vectors in 
K™ [24, Thm 2, 4)]. Since / is continuous on R™, it follows that / is nonnegative 
and convex on R™. If p(x) — it follows that /(y) = on the interval joining x, u. 
Hence <j>{t) = p(x) = 0. Assume that p(x) > 0. Then </>(0) > 0. If 4>{t) is 
constant then <f>(t) = (j>(0) > 0. Assume that <fi(t) is a nonconstant polynomial. Let 
u; > 0, 1 = 1, . . . be a sequence of vectors converging to u. Then <fii(t) :— p(x + tui) — > 
4>{t). Each <j>[ has m real negative zeros. Consider the complex projective space 
CP 71-1 with the homogenous coordinates z = (zi, . . . , z n ) T . The m-roots of 4>i(t) 
correspond to the intersection points of the line Li := {z := sx + tui, s, t € C} with 
the hypcrsurfacc H = {z G CP™ : p(z) = 0} in CP" -1 of degree m. All these points 
are real and are located in the affine part of of the line Li, i.e. s = 1 and t < 0. As 
Li — * L := {z := sx + tu, s,t € C} it follows that L; n H — > £ fl , counting with 
the multiplicities if 0(£) is not constant. Thus all roots of </>(£) are nonnegative. Since 
</>(0) > it follows that all the roots of 4>(t) are negative. Hence 4>(t) > for t > 0. 

Assume that p u ^ 0. Recall that p Ul — > p u . According to part 1 p Ul (y) > p Ui (x) > 
for each pair y > x > 0. Hence p u (y) > Pu(x) > 0. If p u (y) = for some y > 
then for any x > 0, there exists t > such that y > ix > 0. Hence t m_1 p u (x) = 
p u (£x) = p u (x) = => p u = contrary to our assumption. Thus p(y) > 0. As 
Pu; (w + ty) — > V-'(i) : = Pu(w + ty) and ip(t) is a polynomial of degree m—1, it follows 
that ip(t) has m—1 real roots for each w € R™. Thus p u is positive hyperbolic. 

3. Since p(y) > p(x) > for y > x on R"_ it follows that q(yi) > g(x x ) > on R" _1 , 
where y = (y^jO^x = (x^O) 1 ^. Since q is a homogeneous of degree m, from the 
arguments in part 2 it follows that either q{y\) > for each yi > or q = 0. Assume 
that p(y) = q(yi) > 0. Then according to part 1, p(x + ty) = q(xi + tyi) has m real 
roots in t. Assume that r ^ 0. Hence p e „ ^ 0. So p e „ is positive hyperbolic by 2, and 
by previous arguments r is positive hyperbolic on R ra_1 . 

□ 

The following propositions are well known and we bring their proof for completeness. 
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Proposition 2.2 Let p : R" ->Klea positive hyperbolic polynomial of degree m > 1. 
Then the coefficient of each monomial in p is nonnegative. 

Proof. Let p be a positive hyperbolic polynomial of degree m > 1. Let u > 0. Part i 
of the above lemma yields that p u is positive hyperbolic of degree m — 1. 

We prove the proposition by induction on to. Let m = 1. Assume that u > 0. Then 
Pu = Vpu > 0. Hence Vp > and the corollary holds. 

Assume that the result hold for m = I > 1. Let pbea positive hyperbolic polynomial of 
degree to = I + 1. Let = (<5ji , . . . , <5m) T G R™ for i = 1, . . . , n. Let u 3 - > 0,j = 1, . . . , and 
assume that limj—nx, Uj = e^. Hence p Uj is positive hyperbolic of degree I. By the induction 
hypothesis the coefficients of all monomials of p u are nonnegative. Let j — > oo and deduce 
that the coefficients of all monomials of p ei are nonnegative. Hence the coefficients of all 
monomials of p which include the variable Xi are nonnegative. Let i = 1, . . . , n to deduce 
that the coefficients of all monomials of degree one at least are nonnegative. As p(0) > 
we deduce the proposition. □ 



Proposition 2.3 For x = {x u . . . , x n ) T G R" let e x := (e Xl , . . . , e x ") T > 0. Let 
p : M" — > R be a nonzero polynomial with such that the coefficient of each monomial is 
nonnegative. Then logp(e x ) is a convex function on R". Let L C R" be a line in R™. Then 
the restriction o/logp(e x ) to L is either an affine function or a strictly convex function. 

Proof. The convexity of logp(e x ) can be found in [23]. Note that / := p{e XL )\L is of 
the form in Yli=i a i e&i *> where each a* > and b\ > b 2 > ■ ■ ■ > 6fe. If k = 1 then log/(t) is 
logai + bit. Otherwise it straightforward to show that y is increasing on R. □ 

Corollary 2.4 Let p : R™ — > R be a positive hyperbolic polynomial of degree to > 1. 
Then logp(e x ) is a convex function on R". Let L C R ra be a line in R™. Then the restriction 
o/logp(e x ) to L is either an affine function or a strictly convex function. 

Examples of positive hyperbolic polynomials 

1. Let A = (ay )™'J =1 £ ]j mx ™ b e a nonnegative matrix, denoted by A > 0, where each 
row of A is nonzero. Fix a positive integer k g [1,to]. Then 

k 

p fc ,A(x):= X! JJ(Ax) % ,xGR", (2.1) 

l<2i<...ifc<m j-l 

is positive hyperbolic of degree k in n variables. 

2. Let Ai, . . . , A n £ C mxm hermitian, nonnegative definite matrices such that A\ + . . . + 
A n is a positive definite matrix. Let p(x) = det X)"=i x i-^-i- Then p(x) is positive 
hyperbolic. 

Proof. 

1. First note that ^.^(x) > for x > 0. The hypcrbolicity of p m ,A and P\,a is obvious. 
Assume that k G (1,to). Let z = [zi, . . . , z n + m -k) T G R"+ m ~ fc and define P(z) := 

i=i(l^j=i a ij z j Then 

/ m \-i Qrn-kp 

PkM x ) = , a a ((xi, . . . ,a;„,0, . . . ,0)). 

Hence by Lemma 2.1 pfc.A positive hyperbolic. 

2. This is a standard example and the proof is straightforward. 
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Let p(x) : R" — > M be a positive hyperbolic polynomial of degree m > 1. As in [19] 
define 

Cap p := inf p(x) = inf - PW (2.2) 

x>0,Xi...a; n =l x>0 (£1 . . . X n ) » 

It is possible that Cap p = 0. For example let p = x™ 1 ...a;™" where mi,...,m n are 
nonnegative integer whose sum is m and (mi, . . . , m n ) =/= kl. 

Proposition 2.5 Let A £ M. nxn be a doubly stochastic matrix. Let pk,A,k £ [l,n] be 
positive hyperbolic defined as part 1 of the above example. Then Cap pk.A = (jj) ■ Let B £ R n 
be a matrix with positive entries. Then there exists two positive definite diagonal matrices 
Di, D2, unique up to tD\, t~ 1 D2, t > 0, such that A := D1BD2 is a doubly stochastic matrix 
[30]. Let Pn.B be defined as above. Then Cap p n ,B — dct p Do ■ 

Proof. Consider first p n ,A- Since A is row stochastic p n .A(l) = 1- Hence Cap p n ,A < 1- 
Let u = (iti, ■ ■ ■ ,u n ) T be a probability vector. Then for any x = (xi, . . . ,x n ) > 
the generalized arithmetic-geometric inequality states u T x > Yl7=i X T • Use this inequality 
for each (Ax)j. The assumption that A is doubly stochastic yields that p n ,A > Xi . . .x n 
Cap Pn.A > 1- Hence Cap p n .A = 1- 

Let k £ [l,n). Thcnpfc ,a(1) = (")■ Hence Hence Cap Pk,A < (?)■ Apply the arithmetic- 
geometric inequality to %y to deduce that pk,A > (£)PnA- Hence C&ppk,A > (?). 

It is straightforward to show that P ^ B ^ — detTD^D^yl — y~i where y = -D 2 ^ 1 x. Hence 
Capp n , B = dct ^ lD2 - □ 

The following result is taken from [19]. 



Lemma 2.6 Let k > 1 be an integer, u := [u\, . . . , Uk) T > 0, v := (vi,..., Vk) T > 
and define f(t) := Y[i=i( u i^ + K{f) :== m ft>o 27ien /'(0) = A' /or fe = 1 and 

/'(0) > (^-) fc - 1 A^ /or fc > 2. For fc > 2 eguaZ% ZioZds z/ and only if ^ = . . . = ^. 



Proof. We can assume WLOG that /(0) = 1 ; i.e. that f(t) := U.i=i( a i t + 1 ) ^ a ^ = If- 
Using the arithmetic-geometric means inequality we get that 

Kt < f{t) < p(t) = : (1 + ^t) k , t > 0. 

k 

Therefore , by doing basic calculus , 

K<M^- = /'(0)(-A-) fc -\ 
~ t>o t J v Jy k - l' 

which finally gives the desired inequality 

/'(0)>(^) fe - 1 A,fc>2. 

It follows again from arithmetic-geometric means inequality that the equality holds if 
and only if 01 = ^ = . . . = a k = □ 

Definition. Let p : M" ->Rbca positive hyperbolic polynomial of degree m. For each 
integer i £ [0, n] the i-th degree of p is the integer r^ £ [0, m] such that 

Q r ip d r ' p 

— — (xi, . . . ,Xi-i,0, Xi+i, . . . ,x„) ^ 0, and +1 (xi, . . . , x,_i, 0, x i+1 , . . . ,x n ) = 0. 
ox i OXf 

Let deg^ := r.; for i = 1, . . . ,n. 

The following proposition follows straightforward from part 3 of Lemma 2.1. 
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Proposition 2.7 Let p : K™ — > R &e a positive hyperbolic polynomial of degree m. Let 
i G &e an integer. Then 

1. deg;p = p(x) = (p(xi,...,x;_i,0,x; + i,...,x n )). 

i?. For each integer j ^ [0,degjp] ^-f(xi, . . . , St— 1, 0, Xj+i, . . . , x„) is a positive hyperbolic 
polynomial of degree m — j . 

3. For each integer j G [l,n],j =/= i, 

dp 

degj — (xi, . . . ,Xj_i,0,x i+ i, . . . ,x„) < min(degj p, n - 1). 

The following result is crucial for the proof of a generalized Fricdland-Tverbcrg inequality 
and is due essentially to the second author in [19]. 

Lemma 2.8 Let p : W 1 — > M fee o positive hyperbolic polynomial of degree m > 1. 
Assume that Cap p > 0. T/ien deg^ p > 1 /or i = 1, . . . , re. For m = n > 2 

Cap ——(xi, . . . , x,_i,0, x, + i, . ..,£„) > ( — t~ ) deSiP_1 Cap p for i = 1, . . . ,n, 

oxi degj p 

where 0° = 1. 

Proof. It is enough to prove the result for i = re. Suppose to the contrary that p 
does not depend on x n . Then let x(£) = (1, . . . , 1, t) T and £ — > oo in (2.2) to deduce that 
Cap p = contrary to our assumption. 

Assume that m = n > 1. Let fc = deg n p > 1. Let xo := (xi, . . . , x„_i, 0) ,Xi := 
(xi, . . . ,x„_i) T . Proposition 2.7 yields that ff(xi) := |jjr(xo) is a positive hyperbolic func- 
tion in n — 1 variables of degree m — I. Hence ,o(xi) > for xi > 0. Thus for xi > 

k 

p(x + te n ) = k\g(x!)t l + ... = klgfa) JJ(i + Ai(xi)), A;(x) > 0, for i = 1, . . . , k. (2.3) 

i=i 

The second equality follows from part 2 of Lemma 2.1. Assume in addition that x\ . . . x„_i = 
1. Then inf t> o p ( x °+ te ") > Cap p. Apply Lemma 2.6 to the right-hand side of (2.3) to de- 
duce that ^-(xo) > (^-!-) fc_1 Cap p. Since we assumed that x\ . . . x„_i = 1 it follows that 
Cap^(x ) n >(^ i ) fe - 1 Capp. □ 

Remark 2.9 Lemma 2.6 , which is simple but crucial , is a particular case of the 
following general result : 

Let f : [0, oo ) — > J2+ be a nonnegative function differ entiable at zero from the right ; K = 
inft>o ^f-- Ifk>l and fi is concave then /'(0) > (^7— ) fc_1 K . On the other hand if k > 1 
and fi is convex then /'(0) < (^i)*-'" 1 ^. 

3 Friedland-Tverberg inequality 

Theorem 3.1 Let p : R n — > K be positive hyperbolic of degree m G [1, n]. Assume that 
deg^p < r, G [1, m] for i = 1, . . . ,n. Rearrange the sequence r±, . . . ,r n in an increasing 
order 1 < r* < r% < • • • < ?'* • £e£ fc G [1, n] be the smallest integer such that r%. > m — ft. 
Then 
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s ^— ^ (0) - 

l<ii<...<i m <n 1 

f ""v f '"^it-'-^ ![ t r " ' 1 )^-"-Cap » (3.1) 
(n — m ! In — k + 1 ) n rj + n — m 

(Here =1, and the empty product for k = 1 is assumed to be 1.) //Cap > and r{ = m 
for i = 1, . . . , m equality holds if and only if p — C( a:i+ '^ +;r " ) T " /or eac/i C > 0. 

Proof. Suppose that Capp = 0. Then part 3 of Lemma 2.1 yields that the left-hand 
side of (3.1) is nonnegative and the theorem holds in this case. 

Clearly, it is enough to assume the case Cap p = 1. The case m = n is essentially proven 
in [19] and we repeat its proof for the convenience of the reader. Permute the coordinates of 
x\,...,x n such that dcg, n p = min.^^ „] dcg i p < r*. Assume that deg n p = I. Then Lemma 
2.8 yields that r((xi, ■ . ■ , x n -i)) '■= ■§£-{{ x \, ■ ■ ■ s %n-i, 0)) is positive hyperbolic of degree 
n—1 and Cap r > (^p-)' _1 Cap p. Since the sequence (^t^) 1-1 , i = 1, . . . , is decreasing to 
have the lowest possible lower bound we have to assume I = r*. Suppose first that r\ = n. 
Repeating this process n times we get that 



3 — 2 

This inequality to corresponds to the case r* = n for i = 1, ...,n. The equality case is 
discussed in [19]. 

Let m G [l,n - 1]. Put P(x) = pWliElLi^r" 1 - Clearly, P is positive hyperbolic 
of degree n. Since — x i — ( x i ■ ■ ■ x n) ™ f° r each x > 0, it follows that Cap P > Cap p. 
Apply (3.1) to P for m = n to deduce (3.1) in the general case. Since the equality case for 
P holds if and only P = (- X)™=i x i) n ^ follows that the equality in (3.1) holds if and only 

Let A £ M. nxn be a doubly stochastic matrix. Apply this theorem to p m .A defined 
Proposition 2.5 to deduce the Friedland-Tverberg inequality for the sum of all m x m 
permanents of A: 

Corollary 3.2 Let A G R" xn be a doubly stochastic matrix. Then perm m A > {^f ' ^ 
for any m G [2, n]. equality holds if and only A = J n . 

Theorem 3.3 (Gurvits) Let A G M" XTl be a doubly stochastic matrix, such that each 
column contains at most r G [1,ti] nonzero entries. Then 

r\ ,rj-l,(r-i)(n-r) r! , jr . r ( r _i) ,r- -1 . ( r -i)„ 
pcrmA>-(— — ) =_(__) (__) . (3.2) 

Proof. Note that for p(x) = n™=i( J ^ x ) i we have that deg 4 p = r for i = 1, . . . , n. Apply 
(3.1) to this case, i.e. m = n,r* = r, j = 1, . . . , n and k = n — r+l to deduce the theorem. □ 



4 The ALMC and AUMC 

Let G = (V, E) be a general graph with the set of vertices V and edges E. A matching in G 
is a subset MCE such that no two edges in M share a common endpoint. The endpoints 
of the edges in M are said to be covered by M. We can think of each edge e = (u,v) £ M as 
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occupied by a dimer, consisting of two neighboring atoms at u and v forming a bond, and 
of each vertex not covered by M as a monomer, which is an atom not forming any bond. 
For this reason a matching in G is also called a monomer- dimer cover of G. If there are no 
monomers, M is said to be a perfect matching. Note that if a perfect matching exists then 
ffV is even. A matching M with #M = k is called an k-matching. We denote by 4>c{k) be 
the number of fc-matchings in G (in particular 0g(O) = 1), and by $g(x) '■— J2k 4>G(k)x k the 
matching generating polynomial of G. It is known that all the roots of matching polynomial 
are real negative numbers [27]. 

Let G be a bipartite graph G = (V,E), where V = Vi U V2 is the set of vertices of 
G and E is the set of edges that connect vertices in V\ to vertices in V2. Assume that 
#Vi = #V2 = n. We identify V\ and V2 with (n) := {1, . . . , n}, where the vertices in V\ 
and V2 are colored with colors black and white respectively. Then G is represented by — 1 
n x n matrix A(G) = A = (o»j) G {0, l} nXn , where o»j = 1 if and only if the black vertex i is 
connected to the white vertex j. It is convenient to consider multi bipartite graphs. Thus, 
the entries of the representation matrix A(G) = (a^) £ Z" xn are nonnegative integers, 
where a%j is the number of edges from the black vertex i to the white vertex j. 

It is straightforward to show that 

4> G (k) = perm k A(G), k = 0,...,n, where perm A := 1 for any A £ R™ xn . (4.1) 

Let r(n, r) be the set of bipartite r-regular multi-graphs, (each vertex has degree r), with 
n := Let A(n,r) be the set by an n x n nonnegative matrices A with integer entries, 
such that the sum of each row and column is r. Then each G € T(n,r) is represented by 
A G A(n,r) and vice versa. Note for each A £ A(n,r) the matrix - A is doubly stochastic. 
Corollary 3.2 yields: 



Note that the symmetric group S n on n elements, presented as the group of permutation 
matrices II„ C {0, 1}™ X " acts from the left and from the right on A(n, r), i.e. PA(n,r) = 
A(n,r)P for each P £ LT n . These actions are equivalent to the action of S n on V\ and V% 
respectively. 

There is a standard probabilistic model on r(n,r), which assigns a fairly natural prob- 
ability measure v(n,r) on T(n,r) [27]. The measure i/(n,r) is invariant under the action 
of S n on V\ and V2 as explained above. By abuse of the notation we view v(n, r) also a 
probability measure on A(n, r), which is invariant under the left and the right action of n„. 
The following result is proven in [11]: 

Theorem 4.1 Let v(n,r) be the probability measure defined above. Then 




for any G £ F(n,, r). 



(4.2) 



E l y(n^r){4'{G 1 m)) = B 1/( „ iT .)(perm m A) 




m = 0, ... ,n 



(4.3) 



(rn)\ 




-00 



n 



(4.4) 



where 



gh r {p) := \ (plogr-plogp- 2(1 -p)log(l -p) + (r -p)log (l - f)) , 



(4.5) 



The case m = n in (4.3) is given in [27]. 
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Fix as subset J c (n) of cardinality m: f^J = m. For G G T(n, r) let </>(G, J) be all 
TO-matching in G that cover the set J C V 2 . For A = (ay) € M nx ™ and 7 C (n) let A[J|J] 
be the submatrix (a^Oie/.jej- Denote 

perm OT A[(n)|J] = ^ perm m A[J| J]. 

IC(n),#I=m 

Then </>(G, J) = perm m A(G)[(n)\J]. 

Use the invariance of v(n, r) under the action of S n on V% and the fact that there are 
(") distinct subsets J C (n) of cardinality m to obtain: 

Corollary 4.2 Let f(n, r) be the probability measure defined above. Then for any J C 

(n), jfj = m, 

( n )r 2m mHrn - mV 

E„ M ((l>(G,J))=E v(ntr) (perm m A[(n)\J]) = ^ -\. '-, m = 0,...,n. (4.6) 

The following conjecture is stated in [10]. 

Conjecture 4.3 (The Asymptotic Lower Matching Conjecture) 

For r > 2, let G„ = (V n , E n ), n = 1, 2, . . . be a sequence of finite r-regular bipartite graphs 
with #V n — > oo. Let k n G [0, ^|^], n = 1, 2, .. . be a sequence of integers with lim n ^ x ^jj^- = 
p G (0, 1] . TTien 

li m iuf^%^>^( P ). (4.7) 



n — >oo 



For r = 1 this conjecture holds trivially. For r = 2 this conjecture is proved in [11]. The 
inequality (4.2) implies that under the conditions of Conjecture 4.3 the following inequality 
holds, see [13] 

. , log0 G „(fc») 

hmmf — > /ftr(p), (4.8) 

n->oo 

where 

/Ar(p) := i(-plogp-2(l-p)log(l-p)+plogr-p). (4.9) 

As usual, we denote by R[x] the algebra of polynomials in x with real coefficients, by 
G R[x] the zero polynomial, and by M+[x] C R[x] the subalgebra of polynomials with 
non-negative coefficients. We partially order R[x] by writing, for /, g G R[x], g >z f when 
g - / G R+ [x] , and 5 >- / when g - f G K+ [x] \ {0}. Clearly, if 51 ^ /1 ^ and g 2 t /2 ^ 0, 
then >- /1/2 unless gi = fi and .92 = /2- 

Let qK r ^ r denote the union of q complete bipartite graphs K r . r having r vertices of each 
color class. It is straightforward to show that any finite graphs G, G' satisfy 

$GuG'(x) = $ G (x)$ G '(x), (4.10) 

and that 

fc=0 ^ ' 

The following conjecture is stated in [10] 

Conjecture 4.4 (The Upper Matching Conjecture) Let G be a bipartite r-regular 
graph on 2qr vertices where q, r > 2. Then $q ^ ^aK r r > equality holding only if G — qK r r . 

For k = 2 this conjecture is proved in [11]. The above conjecture implies the following 
Asymptotic Upper Matching Conjecture [10]. Denote by K(r) be the countably infinite 
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union of K r ^ r . Let Px( r )(t) and h,K( r )(jp),P G [0,1] be the pressure and the p-matching 
entropy associated and the with K(r) [14]: 



where 



loe:V' n fc'e 2fe * 
W*)= 2, , teK. (4.12) 

h K (rMt)) = P K(r) {t) - tp(t), i G R (4.13) 

^'^wW' f^'XC ' <eK - (4 ' 14) 

Conjecture 4.5 (The Asymptotic Upper Matching Conjecture) 

For r > 2, Zet G„ = (V^, -EVi)> n = 1, 2, . . . &e a sequence of finite r-regular bipartite graphs 
with #V n — > oo. Lei fc„ G [0, ], n = 1, 2, . . . 6e a sequence of integers with limn^^^r- = 
p G (0, 1] . l%en 

,. log <}>G n (kn) . , i , ..... 
hmsup — < h K[r) [jp). (4.15) 

Equality case holds for the sequence qK r .r, 9=1,2,.... 

For r = 2 the AUMC is proven in [11]. For p = 1 and any reff the AUMC follows from 
the proof of Mine conjecture by Bregman [4]. Some computations performed in [10] support 
the ALMC and AUMC. 

The following plots illustrating the Asymptotic Matching Conjectures for r = 4, 6. Let 
C n a cycle on n points, and let T n ^ = (V n , E n ) := C„ x . . . x C n , n = 3, ... be a sequence 

d 

of d dimensional torii. Note that each T n ^d is 2d regular graph. It is a classical result that 
the following limit exists for any p G [0, 1]: 

llm = h d(P), P G [0, 1]. (4.16) 

hd(p) is the d-dimensional monomer-dimer entropy of dimcr density p G [0, 1] in the lattice Z d 
[21] and [13]. In this case we use the notation ha '■= rnaXp 6 [o,i] hd(p), (the quantities /id and 
hd := ftd(l) are called the d-monomer-dimer entropy and the 2-dimer entropy, respectively, 
in [13]). Figure 1 shows various bounds and values for the monomer-dimer entropy /12O?) of 
dimer density p G [0, 1] in the 4-regular 2-dimensional grid. FT is the Friedland-Tverberg 
lower bound fh±(jp) of (4.9), h2 is the true monomer-dimer entropy equal to max pe [ .i] h^ijp) 
(it is known to a precision much greater than the picture resolution). The crosses marked B 
are Baxter's computed values [1]. ALMC is the function gh±(jp) of (4.5), conjectured to be a 
lower bound in the Asymptotic Lower Matching Conjecture. AUMC is the monomer-dimer 
entropy hx(p) of dimer density p in a countable union of ^4,4, given by (4.12)-(4.14) and 
conjectured to be an upper bound by the Asymptotic Upper Matching Conjecture. Notice 
that AUMC goes a little over h2: a countable union of ^4,4 has a higher monomer-dimer 
entropy than an infinite planar grid. 

Figure 2 shows similarly various bounds and values for the monomer-dimer entropy /13 (p) 
of dimcr density p G [0, 1] in the 6-regular 3-dimensional grid. FT is the Friedland-Tverberg 
lower bound fhs(p) of (4.9), h3Low and h3High are the best known lower and upper bounds 
for the true monomer-dimer entropy equal to max pe [ .i] h^{p). ALMC is the function gh§{p) 
of (4.5), conjectured to be a lower bound in the Asymptotic Lower Matching Conjecture. 
AUMC is the monomer-dimer entropy Hk{p) of dimcr density p in a countable union of 
-^6,6, given by (4.12)-(4.14) and conjectured to be an upper bound by the Asymptotic 
Upper Matching Conjecture. Notice that AUMC goes a little over h3High: a countable 
union of Kq^ has a higher monomer-dimer entropy than an infinite cubic grid. 
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P 

Figure 1: Monomcr-dimcr tiling of the 2-dimcnsional grid: entropy as a function of dimcr 
density. FT is the Friedland-Tverberg lower bound, h2 is the true monomer-dimer entropy. 
B are Baxter's computed values. ALMC is the Asymptotic Lower Matching Conjecture. 
AUMC is the entropy of a countable union of i^4,4, conjectured to be an upper bound by 
the Asymptotic Upper Matching Conjecture. 



5 A proof of some case of the ALMC 

In this section we prove the following case of ALMC: 

Theorem 5.1 Let r > 3 be an integer. Then the asymptotic lower matching conjecture 
(4.8) holds for Ps = ^r s ,s = 0,1,2,... . 

The proof of this theorem follows from the following results. 

Theorem 5.2 Let p : R n — > R be positive hyperbolic of degree m G [1, n). Assume that 
deg^ < ri G [l,m] fori = l,...,n. Rearrange the sequence ri,...,r„ in an increasing 
order 1 < r\ < r^ < .. . < r*. Let s £ N. Let k € [l,n] be the smallest integer such that 
r* k + s > n — k. Then 

s dXn .X w ^ 

l<ii<...<i m <n 1 

- — ( T fc t 1) L i n(^ +g ' 1 ) r - +S - 1 ^P P . (5.1) 

s n -' m (n - m)!((s - l)n + m)\ (n - k + 11 v r* + s ' y ' 

Proof. Let q : R™ — > R be positive hyperbolic of degree n — m with deg 4 g < s for 
i = 1, . . . , n and Cap q = 1. Then / = : R" — > R is positive hyperbolic of degree n with 
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0.2 0.4 0.6 0.8 1 

P 

Figure 2: Monomer- dimer tiling of the 3-dimcnsional grid: entropy as a function of dimcr 
density. FT is the Fricdland-Tverberg lower bound, h3Low and h3High are the known 
bounds for the monomer-dimer entropy. ALMC is the Asymptotic Lower Matching Conjec- 
ture. AUMC is the entropy of a countable union of Kq.q, conjectured to be an upper bound 
by the Asymptotic Upper Matching Conjecture. 



Cap / > Cap p and deg^ / < rj + s for i = 1, . . . , n. Apply Theorem 3.1 to / to deduce 

am an—m 

y d_p d q > 

\ j=l J 

where l<i[ < ... < i' n _ m < n and i m , i[, i' n _ m } = (n). 

Let A £ A(n, s) and choose q = ( ™ ) P n _ m Ia( x ) as m (2-1). Note 



. . . dxj 



■(°) = 7 n r^ perm n _ ro i[{n}l/] 

vn— m/ 



Now take the expected value of the left-hand side of the inequalities (5.2) corresponding to 
all A G A(n,s). Use Corollary 4.2 to deduce that the coefficient of each Qx d ^ x (0) is 

~ m (n— m)! ((s— l)n+m) ! |— j 



(sn)! 
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Corollary 5.3 Let p : R n — > R be positive hyperbolic of degree m € [1, n). Assume that 
degiP <r£ [1, m] for i = 1, . . . , n. Let s G N and k — n — r — s + 1 > 1. TTien 



l<ii <. . .<z m <n 

(sn)! (r + s)! /r + s — 1 



y — (°) ^ 



l (n - m)!((s - l)n + m)! (r + s) r+s 



)(r+ S -l)(„-r- S ) Cap ^ (53) 



Theorem 5.4 Let B G R" xn &e a doubly stochastic matrix with at most r nonzero 
entries in each column. Let s G N and k = n — r — s + 1 > 1. Then for each m G 

perm m B > , - ^ (r^ r_ + s-l (r+s _ 1)(n _ r _ s , 

~ s n - m (n-m)!((s- l)n + m)! (r + s) r+sV r + s ; v ; 

Proof. Let p = p m ,s(x) as defined by (2.1). Then (5.4) follow from Corollary 5.3. □ 

Let G„ G r(n, r). Then G„ is represented by its incidence matrix A n G A(n, r). Let 
_B„ := ^A„. Then i3 n is a doubly stochastic matrix where each row and column of B n has 
at most r positive entries. Clearly, the ALMC conjecture follows from the following stronger 
conjecture: 

Conjecture 5.5 (The Asymptotic Lower r-Permanent Conjecture ) 

For r > 2, let B n , n = 1, 2, . . . be a sequence of nx n doubly stochastic matrices, where each 
column of each B n has at most r-nonzero entries. Let k n € [0, n], n = 1, 2, . . . be a sequence 
of integers with lim n ^oo^ = p G (0,1]. Th 



en 

, log perm fc B n p 
hmmf — 2 > gh r (p) - - logr. (5.5) 

n— too In I 

Theorem 5.1 follows from the following result: 

Theorem 5.6 Let r > 3, s > 1 be integers. Let B n , n = 1,2, ... be a sequence of n x n 
doubly stochastic matrices, where each column of each B n has at most r-nonzero entries. 
Let k n G [0, n], n = 1, 2, . . . be a sequence of integers with lim n ^oo = P £ (0, 1]. Then 

Uminf lQ gP erm k n B n ^ 1 (_ plogp _ 2 (1 - p) log(l - p)) + (5.6) 

ri— >oo in A 

\ ((r + s-l) log(l - -L-) -(s-l+p) log(l - 
2 \ r + s s 

Moreover, the Asymptotic Lower r -Permanent Conjecture 5.5 holds for p s = ;q-;,s = 
0,1,2,.... 

Proof of Theorem 5.6. Apply the inequality (5.4) to B n for m = k n . Take the 
logarithm of the both sides of this inequality and let n — > oo. A straightforward calculation 
for the right-hand side, using the Stirling's formula, yields the inequality (5.6). Assume 
that p = p s = Then ~ Ps = = £f. Then the right-and side of (5.6) is equal 

to gh r (p s ) — ^logr. Hence the asymptotic lower r-permanent conjecture 5.5 holds for 
p a = ^j,* = l,2,.... 

We now discuss the case s = 0, i.e. p = p a = 1. Let B = (by)"j = i be any n x n 
nonnegative matrix. Denote by G(B) = (V,E) the bipartite graph induced by B, i.e. the 
edge (i,j) is in E, if and only if bij > 0. Then B induces the weighted graph on G, where the 
weight of the edge isbij. Let pb(x) = x n + ^" i=1 (— l) m perm m (B). ps(x) is called the 
matching polynomial of the weighted graph G. Heilmann and Licb showed in [22] that Pb{%) 
has nonnegative roots. (See also [27].) Hence the arithmetic-geometric inequality for the 
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elementary symmetric polynomials of the nonnegative roots of pb{x) yields the inequality 
perm m B > ( "J (perm B) 1 ^ . (See [34] for the case of m-matchings in bipartite graphs.) 

Use Theorem 3.3 to deduce that perm B n > ^(^ziY^ ^i 1 -^-)^ ■ Apply the above 
two inequalities for the sequence B n and m = k n for n — 1, 2, ... to deduce the case p — 1. □ 

Let d = 6 and p 3 := f = §. Then Theorem 5.1 yields that /i 3 (§) > .7845241927, which 
implies that h 3 = max pg [ 0i i] h 3 (p) > .7845241927. This improves the lower bound implied 
by (4.8) h 3 > .7652789557 [13]. The computations in [13] yield that h 3 < .7862023450. 
Thus h 3 S [.7845,.7863]. 



6 Matching in general graphs - hafnians 

Let G = (V, E) be a graph on the set of vertices V and the set of edges E. Assume that 
jfV = N. Then G is represented by a symmetric 0—1 matrix B = B(G) with diagonal. 
If G has a perfect matching then N = 2n is even. If G is bipartite and V = V\ U V-z , where 
V\ = {1, . . . , n}, V% = {n + 1, . . . , 2n}, we deduce that 

B =0 o)' (61) 

where A is the representation matrix of the bipartite graph G. As explained above, the 
number of m-matching in the bipartite graph G is perm m A. 

In this section we discuss the m-matching of a general graph G, and the related function 
haf m B which counts the number of m-matching in G. Let a, P C {1, 2, . . . , N} be two 
nonempty sets of cardinality #a = = j, respectively. We then arrange the 

elements of a = {ai, . . . , o^} and (3 = {/3i, . . . , /3j} in an increasing order: 1 < a,\ < . . . < 
on < N, 1 < /3i < . . . < fa < N. For B = (b st ) e C NxN we denote by B[a\0\ € C lX] 
the submatrix (6 Qsft )^ =1 . Denote by Sj(R) D S;(M+), S ij0 (R) D Si, (K+) the space of real 
valued I x / symmetric matrices, the cone I x / symmetric matrices with nonnegative entries, 
the subspace of real valued I x / symmetric matrices with zero diagonal, and the subcone 
of Z x / symmetric matrices with zero diagonal and nonnegative entries respectively. Let 
B G Sjv(R) and an integer m £ [1, -y]. Then the m — th hafnian of B is defined as 

haf TO B = 2" m permB[a,P]. (6.2) 

a,/3c{l,...,iV},#Qi=#/3=)7i,Qin/3=0 

That is if ■ ■ ■ , {imjjm) is an m matching of a complete graph Kn on vertices, 

then the product 6, 131 . . . bi m j m appears exactly once in haf m B. Since bi t j t = bj l i l there are 
2 m choices of a and /? for which this product appear, we need to use the factor 2~ m in the 
above definition of haf m B. If B = B(G) then haf m B gives the number of m-matching 
in G. Note that from the definition of perm m B it follows that haf m B < 2 _m perm m B. 
Equivalently, it is straightforward to show: 

haf m B = (2 m m!)- 1 ^ (x T i?x) m , BeS N (R). (6.3) 

i<i 1 <...< l2m <w ° Xil ■■■ ax ^ 

Unfortunately, the quadratic polynomial x T 5x is not always positive hyperbolic. Note that 
haf m B does not depend on the value of the diagonal entries B. Let B^ be the matrix 
obtained from by replacing the diagonal entries of B by zero elements. Then haf m B = 
haf m £(°). 

For B £ S;(M) we denote by \i(B) > ... > M(B) the I eigenvalues of B counted 
with their multiplicities and arranged in the decreasing order. As usual for B,C G S;(M) 
we denote C ^ B if C — B is a nonnegative definite matrix. The maxmin, or minmax 
charaterization of Xk{B) yields that if C >z B then Afc(C) > Xk(B) for k = 1, ...,Z. In 
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particular, if B £ S;(R) has nonnegative diagonal entries then B y B^ and \k(B) > 
X k (B^)fovk = l,...,l. 

The following result is well known and we bring its proof for completeness. 

Lemma 6.1 Let B £ S„(R) and n > 2. Then x T Bx is positive hyperbolic if and only 
ifO^B £ S„(R+) and > X 2 (B). 

Proof. Assume that x T £>x is positive hyperbolic. Then Proposition 2.2 yields that 
B £ S„(R+). Since x T Bx > for each x>05/0. Hence Xx(B) > 0. 
Observe next that the positive hyperbolicity of x T i?x is equivalent to 

(x T 5y) 2 > (x T Bx)(y T By), for any x > 0,y £ R™. (6.4) 

Clearly, the above condition holds for any x > 0. The Perron-Frobenius theorem yields 
that there exists 0/xe M™ such that Bx = Xi(B)x. Let ^ y £ R™ and assume that 
y T x = 0. Then (6.4) yields that y T By < 0. Hence X 2 (B) < 0. 

Vice versa suppose that ^ B £ S n (R+) and X 2 (B) < 0. Recall that there ex- 
ists a permutation matrix P £ {0,l} ,lXTl such that P 1 BP is a block diagonal matrix 
diag(£>i, . . . , Bk), where Bi £ S ni (R+) is irreducible. So X\(Bi) > unless m = 1 and 
Bi = 0. Hence our assumptions yield that we may assume that / Bi 6 S ni (R + ) is 
irreducible and Bi = for i = 2, . . . , k. Thus it is enough to show that x T Bx is positive 
hyperbolic for an irreducible B £ S„(R + ), where n > 2. Clearly x T Bx > for x > 0. Thus 
it is left to show that (6.4) holds. Assume first that 6x = Ai (B)x, x > 0. Then (6.4) follows 
straightforward. Suppose x > is any vector. Then there exists a unique diagonal matrix 
D, with positive diagonal entries such that Bx = D~ 2 x. That is DBD(D~ 1 x) = D~ 1 x. Re- 
placing B, x, y by DBD, D~ 1 x, Z? _1 y we deduce the inequality second inequality of (6.4). □ 

Definition 6.2 Let G = (V, E) be a graph on the set of vertices V . Then for k > 2 G 
is called k-partite, if the following condition holds. There exists a decomposition of V to a 
disjoint union of k nonempty sets V\ , . . . , such that each edge e £ E connects Vi to Vj 
for some i ^ j . G is called a complete k-bipartite, if there exists a decomposition of V to a 
disjoint union of k nonempty sets V\ , . . . , Vu such that E consists of all edges from Vi to Vj 
for all 1 < i < j < k. 

Note that a complete graph G with n > 2 vertices is complete n- partite. 
We need the following elementary result whose proof is straightforward. 

Proposition 6.3 Let n > 3 and F = [fij] £ S n (K) with V s on the main diagonal. 
Suppose that and all subsets of cardinality three a = {l<i<j<k<n} the principal 
submatrix F[a,a\ is nonnegative definite . If fij = 1, fjk = 1 then fik = 1. 

Lemma 6.4 Let B = [bij] £ S n (R+), i.e. B is a real n x n symmetric matrix with 
nonnegative entries. Denote by G(B) = ((n),E) the graph, (with no self-loops), induced 
by B, i.e. £ E if and only if i ^ j and bij > 0. Assume that B is irreducible and 

X2(B) < 0. Then G(B) is a complete k-partite graph for some k £ [2,n]. 

Proof. Since £?' ' is irreducible and B y B(°\ it is enough to prove the lemma 
in the case where all the diagonal entries of B are equal to zero, i.e. B = B(°>. Let 
x = (xi, . . . , x n ) T be the unique positive eigenvector of B, corresponding to the maximal 
eigenvalue X\(B), satisfying the condition x T x = Xi(B). Then A = xx T — B and A y 0. 
Let D = diag(xi, ...,x n ). Then the zero pattern of the matrix C = D~ 1 BD~ 1 is the same 
as of the matrix B. Let F = [fij] := 11 T — C , where 1 := (1, . . . , 1) T is the vector of all 
ones. Then F — [fij] > 0. Notice that fa = 1,1 < i < n, and for i ^ j fij = 1 if and 
only bij = 0, (i.e. the vertices i,j are not connected in the graph G(B)). As any principal 
submatrix of F is nonnegative definite, it follows from Proposition 6.3 that for any triplet 
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l<i<j<k<n such that if = 1, fjk = 1 then also fik = 1. In other words, the 
relation "i ~ j" <^=4> "fry = 0" is an equivalence relation. Therefore the graph G(B) is 
complete fc-partite, where each equivalence class of vertices corresponds to some class Vi in 
the fc-partite graph. 

Theorem 6.5 Let G be a graph on n > 1 vertices, and denote by A(G) the incidence 
matrix ofG. Then x T ^4(G)x is positive hyperbolic if and only if G is a union of a complete 
k(> 2) -partite graph on at least two vertices and of isolated vertices. 

Proof. We first show that for a complete /c-partite graph G on at least n > 2 vertices 
"^(G) < 0, which is equivalent to the positive hyperbolicity of x T /l(G)x in view of Lemma 
6.1. Let J n be a symmetric matrix whose all entries are equal to 1. Then J is rank one 
matrix with Ai(J n ) = n and Xi(J n ) = for i = 2,...,n. Let k € [1,^], ni € N,i = 
1, . . . , k, 1 < nk < ... < m and ni + . . . + nk = n. Consider the block diagonal matrix 
J(ni, . . . ,rik) '■= diag(J ni , . . . , J nk )- Clearly J(ni, . . . ,rik) is a nonnegative definite matrix. 
It is straightforward to see that renaming the vertices of G, we will obtain that A(G) = 
Jn — J(fii, . . . , rik) for some unique n\ > . . . > 1. Then minimax characterization of 
X 2 (A(G)) yields that X 2 (A(G)) < A 2 (J„) = 0. Hence x T A(G)x is positive hyperbolic. 

Assume now that x T A(G)x is positive hyperbolic. Therefore G must have at least one 
edge and X2(A(G)) < 0. Hence G has one connected component containing at least two 
vertices and a union of isolated vertices. Without loss of generality we assume that G is 
connected. Then A(G) satisfies the assumptions of Lemma 6.4 and G is fc-partite. □ 



Remark 6.6 Let A be a symmetric nxn matrix with nonnegative entries. It is straight- 
forward to show that the polynomial x T Ax is positive hyperbolic if and only if the function 
V x T Ax is concave on the positive orthant i?" . If A is real, symmetric, nonnegative definite 
then V x T Ax is convex on R n . In view of Remark 2.9 , it is natural to conjecture that if 
A is 2n x 2n real, symmetric, nonnegative definite then the reverse van der Waerden bound 
holds : 

2™n!haf n B < y^Cap (p), p(x) := (x T Ax)". 



7 Algorithmic applications 

One of the main purposes of this paper is to construct a generating homogeneous polynomial 
p(x) of degree n, for some quantity of interest Q as hafnian: haf A, sums of subhafnians: 
haf m A, permanent: perm A, sum of subpermanents: perm m A, such that Q = g x 9 P g x (0) ■ 
If such polynomial is positive hyperbolic then we can apply the results from [19] and the 
results of the previous sections of this paper to get a lower bound on Q, and even to get 
deterministic polynomial-time algorithms to approximate Q within simply exponential factor 
as in [19]. In the general, (not positive hyperbolic case), we can use this representation to 
obtain exact algorithm to compute Q in 2 n poly(n) number of arithmetic operations provided 
the the generating polynomial p can be itself evaluated in poly(n) number of arithmetic 
operations. We present below some examples of this approach. 



7.1 Formula for V 1<5 . , , 5 g% (0) 

Our exact algorithms are based on the following elementary identity 7.2. 

Let p(x) be a polynomial of degree m in n > m variables , p(0) = 0. Define 

Sj= ^2 P(h,-,b m ). (7.1) 

!>je{0,l},£ l<j<m=i 
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Let d„ = (rf n> i, rf ra ,m-i) be the unique solution of the system of linear equations 
d n A — (— 1, — 1) , where the lower triangular matrix A = [o»j] is defined 

as follows: 

TL — j\ 

I h" i > j and = otherwise. 
Then the following equality holds : 

£ a,, a,. (°)=kw) + E ( 7 - 2 ) 

l<ii<...<i m <« 11 lm l<j<m-l 

Notice that this formula requires X]o<i<n-i (j) evaluations of the polynomial p. 
The formula 7.2 follows from the following obvious identities : 

Si= ^2 ( • . ] Cj, 1 < i < m, 
i<j<i ^ •* ' 

where Cj is the sum of the coefficients of all monomials in p involving exactly j variables. 

The formula 7.2 is, in a sense, optimal for n = m, i.e., if for all homogeneous polynomials 
p(x) of degree n 

d n p 

(x) = ^2 a,ip(zi), a i eC,z. l eC n , 



dx 1 ...dx n 

Ki<k 



thenfc>(p«^ [18]. 



7.2 Ryser' like formulas for sums of subhafnians and subperma- 
nents 

1. Let B G Sat(K),x := (xi, . . . ,x N ) T G C w and m G [l,n] n N. Define, as in the proof 
of Theorem 3.1, the polynomial : 

P m (x) = -^(x T Bx)'» 
2 m m! 

As haf rn B = (2 m m!)" 1 Ei<i 1 <...<i 2ra <jv 5^3^— ( xTfix ) m . the application of 7.2 
gives the following Ryscr-like formula for haf m B: 

haf m B = P m {l,. ..,!)+ Yj s i d n,ii ( 7 - 3 ) 

l<j<ra-l 



where Sj are defined by (7.1) forp = P m . The formula (7.3) provides SB(N 7 m)(0(N 2 )- 
0(log(m))) , algorithm to compute haf m B, where SB(N,m) := J2o<j<2m~i (j) • 

2. Forx:= (xi, x n ) T G C" let 



^m(^) ^ ^ ^ii_"'^i m ' 

l<i 1 <...<i m <n 

be the m — th symmetric function of x. Let A be n x n complex matrix and define 
p m (x) := S' m (ylx). Then p m (x) can be evaluated in Oln 2 ) arithmetic operations and 

perm m A= £ — — p m (x). 

l<ii<...<i m <n 

Which gives the following algorithm to evaluate perm m A 

perm m A=p m (l,...,l) + ^ Sjd n ,j, (7.4) 

l<j<m-l 
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where Si are defined by (7.1) for p = p m . The formula 7.4 provides SB(N, m)(0(N 2 )) 
algorithm to compute perm m A. 

Notice that the naive algorithm , i.e. computing and adding all mxm subpcrmancnts, 
requires ((^)) 2 2 m O(m) arithmetic operations . 

7.3 Positive hyperbolic polynomials and convex relaxations 

In this section we always assume that x = (x±, . . . ,x n ) T ,y = (yi, . . . ,y n ) T G C", 1 = 
(1, . . . , 1) T . Suppose that a positive hyperbolic polynomial 

P( x ) = X! a (n,...,r„) n ai[Sm>n 

£l<i<n r ' = " l<i<m 

has nonnegative integer coefficients and is given as an oracle. I.e. we don't have a list 
coefficients , but can evaluate p(x) on rational inputs . The number logp(l) measures the 
complexity of the polynomial p . 

A deterministic polynomial-time oracle algorithm is any algorithm which evaluates the 
given polynomial p(x) at a number of rational vectors qW = (q[ l \ q„ ) which is polyno- 
mial in n and logp(l); these rational vectors q^' are required to have bit-wise complexity 
which is polynomial in n and logp(l); and the number of additional auxiliary arithmetic 
computations is also polynomial in n and logp(l). 

If the number of oracle calls, (evaluations of the given polynomial p(x)), the number of 
additional auxiliary arithmetic computations and bit-wise complexity of the rational input 
vectors q^- 1 are all polynomial in n, (no dependence on logp(l)) then such algorithm is 
called deterministic strongly polynomial-time oracle algorithm. 
The following theorem was proved in [19]. 

Theorem 7.1 There exists a deterministic polynomial-time oracle algorithm, which 
computes for given as an oracle indecomposable positive hyperbolic polynomial p(x) a number 
F(p), satisfying the inequality 

^ -(0)<F (P )<2^^V(0) 



dx\. ..dx n n\ dx\...dx n 

Our goal in this paper is to extend Theorem 7.1 to approximate 
J2i<i 1< ...<i m <n (°) f° r a S iven positive hyperbolic polynomial g(x) of degree m . 

The algorithm behind Theorem 7.1 is based on two observations . First, 

logCapp= inf logp(e y ), (7.5) 

Ei<i<„vi=° 

where we used the notation of Proposition 2.3. If the coefficients of the polynomial p are 
nonnegative then the functional logp(e y ) is convex , the indecomposabilty of the polynomial 
p is exactly uniqueness and existence of the minimum in 7.5 . 
Second point is the inequality : 

Cap p— < P — (0) < Cap p (7.6) 

n" ax\...ax n 

If the positive hyperbolic polynomial p is not indecomposable , wc need first to check if 
Cap p > . If Cap p = then also 9x ^ p dx (0) = . In the case Cap p > we " slightly" 
perturb the polynomial p to get the indecomposability. 

Theorem 3.1 in this paper provides an analogue of the left inequality in (7.6) for positive 
hyperbolic polynomial g(x) of degree m < n . The problem is that in this case we don't 
have the right inequality. I.e. it is possible that ^2 1<i±< <j TO < n g x . ~ g x . — 9(0) > but 
Cap (q) — . This problem can be easily overcome by the following equivalent reformulation 
of Theorem 3.1: 
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Theorem 7.2 Consider a positive hyperbolic polynomial g(x) of degree m < n . Define 

£ x% 

a positive hyperbolic polynomial p(x) = q(x)( — — l) n_m anrf 

(n-m)! \- d m g , n >, 

l<ii<...<i m <n 1 m 

TTien i/ie following inequality holds : 

—Cap p < D m (q) = g x P gx (0) < Cap p. (7.7) 

Corollary 7.3 Let q(x) be a positive hyperbolic polynomial of degree m < n given as 
an oracle. Then there exists a strongly polynomial-time, (in n), oracle algorithm to check if 

El<. 1 <...<v<n dx;, ...ax,- (0) > ^- 



Proof. It follows from Theorem 7.2 that Y J - I<4 , ^ <-„ a (0) > if and 

only if g 9 ^ (0) > , where p(x) = ff(x) ( — 1 ~^ t ~" )"~ m is positive hyperbolic polynomial 
of degree n . Notice the polynomial q is easy to evaluate given an oracle evaluating the poly- 
nomial p . Let q(x) = ^ r i =n a ri,-,r n x i 1 ■■■ x n > > tne support Supp(q) = {(ri,...,r„) : 
Qn,...,r„ 7^ , the Newton polytope is the convex hull of the support CO(Supp(q)) . It 
was proved in [18] that an integer vector (n, ...,r„) G Supp(q) if and only if (ri, ...,r„) £ 
CO(Supp(q)) . Corollary 4.3 in [19] provides a strongly polynomial (in n) oracle algorithm 
for the membership problem " X £ CO(Supp(q))7" for positive hyperbolic polynomial q(x) 
of degree n . 

□ 



Corollary 7.4 Let g(x) 6e a positive hyperbolic polynomial of degree m < n given as an 
oracle. Then 

1. There exists a strongly polynomial-time, (in n), oracle algorithm to check if 
D m {q) := E,< il< ... <im < n axS-kJ ) > 0. 

2. There exists a deterministic polynomial-time oracle algorithm which computes a num- 
ber F m (q) satisfying the inequality 1 < yjl < 2~/(n,m) < 2-^-. 

Corollary 7.5 Let A be n x n matrix with nonnegative entries. Then there exists a 
deterministic polynomial-time algorithm which computes a number P m (A) satisfying the 
inequality 1 < ~ a ^~4 < 27(71, m) < 2^-. 

It s very possible that there exists a deterministic polynomial-time algorithm (in n) which 
approximates perm m A within multiplicative factor e m . 

7.4 A conjecture 

Let the assumptions of Corollary 7.4 holds. To approximate D m (q) we used the identity 

<:)" 



dx\...dx n 



g(0)r(0) = D q {m), (7.; 



for r(x) = ((n — m)\) Q^ILi x i)" m • ^ is natural to ask we can improve our estimates if 
we choose a different positive hyperbolic r(x) such that (7.8) holds. 
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Definition 7.6 Denote the set of positive hyperbolic polynomials of degree m and in n 
variables as PHYP(n,m); and by UPHYP(n,m) C PHY P(n,m) the subset of polynomi- 
als satisfying dx 9 dx . r(0) = 1 for all 1 < i\ < . . . < i n -m < n - Define 

i \ . , Cap qr 
7(77,777) = mi sup ~^ . 

reUPHYP(n,n-m) qeP HYP(n,m) dx *.. dXn <l{Q, ■■■0)r(0, 0) 

Assume that r E UPHYP(n,n — m). Since r is positive hyperbolic, all its monomials 
have nonnegative coefficients. Hence r(x) > 5„.„- m (x) for any x > 0. In particular 
Cap r > Cap S n , n -m = (") ■ By choosing in inf sup definition 7(71, m) q= C ri+ '„' +3: " ) m wc 
deduce straightforward that 7(71, m) > ^y. 

Conjecture 7.7 

Cap Q'S'n.n-m n m 



7(n, m) = sup 



q£PHY P(n,m) Q Xl ^g Xi 

Note that the hyperbolic van der Waerden inequality [19] implies that 7(71,771) < ^-r. It 
follows from Theorem 3.1 that for each q e PHYP(n,m) such that Cap q = q(l) we have 
the inequality 

Cap qS ntn — m 7i m 



dx\ . ..dx. 



-g(0)5 n ,„_ m (0) to! 



Remark 7.8 M^e presented in this section one "natural" generating polynomial for the 
hafnian, and described all symmetric boolean matrices such that this polynomial is positive 
hyperbolic. It is an interesting open problem whether there exists a generating positive hyper- 
bolic polynomial for the hafnians of boolean matrices which can can evaluated in polynomial 
time. If the answer to this problem is negative it can explain why approximating the hafnian 
(number of perfect matchings in general graphs) is "harder" than the same problem for the 
permanent. It is easy to prove that computing the hafnian of integer symmetric 2n x 2n ma- 
trices with nonnegative entries and the signature (+,—,...,—) if #P- complete. The results 
in this paper allow to use Sinkhorn's scaling to approximate the hafnian within multiplicative 
factor e 2n in this "hyperbolic" case. 
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